

i 

N76-15031 


THE PLANAR DYNAMICS OF AIRSHIPS 


Frank J. Regan* 

L 

i 


ABSTRACT ; This paper will consider the forces and moments 
acting upon a LTA vehicle in order to develop parameters 
describing planar motion. Similar expressions for HTA 
vehicles will be given to emphasize the greater complexity 

! of aerodyn«unlc effects when buoyancy effects cannot be 

neglected. A brief summary is also given of the use of 
virtual mass coefficients to calculate loads on airships. 



q 


a 

a 


q 

q 


a 



SYMBOLS 

Drag coefficient 

Pitching moment coefficient/ Ny/QSi 
d Cm/a (qt/2V) 

3 Cm/3 
3 CSV^3a 

3 Cm/3 (4i/2V) 

Normal force coefficient/ Fz/QS 
3 Cz/3 (qt/2V) 

3 Cz/3 (4t^/2V^) 

3 Cz/3a 

3 Cz/3 (ot/2V) 

Drag force 

Gravitational acceleration 
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Transverse moment 
Transverse radius 
Body length 
Reference length, 

Pitching moment 
Body mass 
Dynamic pressure, 

2/3 

Reference area, V ' 

Airship density to medium density, ^b/p 

Airship speed 

Airship volume 

Inertial axes 

’Body axes 

Normal force 

Angle of attack 

Angle of pitch 

Velocity potential 

Density 


of inertia ^ 

of gyration, yiy^I^ 

body length 


1/2 pyo** 


INTRODUCTION 

In studies of the dynamics of Heavier Than Air (HTA) vehicles, effects 
due to buoyancy are almost invariably neglected. Sustaining force is 
the result of relative motion existing between the HTA vehicle (or at 
least some portion of the vehicle) and the surrounding air mass. In 
short, the lift or sustaining force associated with HTA craft is 
entirely dynamic. 

A somewhat reverse situation exists in the case of Lighter Than Air 
(LTA) craft. The principal sustaining force comes from buoyancy, with 
perhaps a small additional force (about 10 percent) available under 
lome conditions from dynamic lift. To put the con^arison between LTA 
and HTA craft on at least a semiguantitative basis, it is convenient 
to define a relative density parameter, s, as 

<i) 

It may be seen that s is of 0 (1) for a LTA vehicle, while for a HTA 
s is no less than 0 (10'*^) and for swat cases 0 (10^*). 

In addition to buoyancy playing an essential role in LTA dynamics, 
there are in addition dynamic effects which for convenience might be 
lumped in the terms virtual mass. Such dynamic effects are ti.ken to 
mean forces and moments arising from (and hopefully linear with) 
angular rate or linear acceleration. These virtual mass effects are 
essentially reactive forces and moments caused by imparting an angular 
velocity and a linear and angular acceleration to the surrounding air. 
Like buoyancy these virtual effects are usually neglected fox HTA 
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craft; for LTA vehicles/ however/ such effects form an essential part 
of the loads acting on the craft. Thus such effects enter prominently 
into any considerations of stability. 

No originality is claimed in the following development of either the 
mathematical model of planar dynamics or the subsequent load calcula- 
tion methods. The equations of planar motion originated with ballisticians 
such as Murphy . However/ because of the negligible effect of 
buoyancy/ great simplifications are possible in the aeroballistic 
formulation. As will be shown/ the airship equations are far more 
complex. The load calculation techniques follow from Bryson '2) origi- 
nally and have been presented by Nielson Again these methods are 
applied to LTA vehicles rather than the HTA missiles which were the 
original motivation for Bryson's work. 

DYNAMICS OF PLANAR MOTION 

Consider an airship luidergoing planar motion as illustrated in 
(1) below 






FIG. 1 FORCES AND MOMENTS ACTING ON AIRSHIP 


The axes Xe/Ye/Ze are the inertial axes, while X/Y,Z are body-fixed 
axes. The equations of planar motion are the forces along axes XS/Ze 
ano the moment about axis Ye. Note that because of the definition o£ 
planarity axis Ye is identical to axis Y. 

The moment and two force equations may be written as 

mVCos<e-«c)- <2a> 






F*Ccse- (\- - 5) (2b) 

I 0 - My ( 2 c) 

where e is the relative density parameter of equation (1) . In addition 
to the three load equations above » figure (1) also provides the follow- 
ing kinematic relationship: 

* -T (3a) 

which gives upon differentiation 

VSinW-«<)-VC.stt-‘)te-«Q (3b) 


Under the assumptioa that the X axis does not greatly vary from the 
horizontal, Xe, it is possible to restrict e and a to small angles. 
Subject to such small angle restrictions equations (2) and (3) become: 
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(4c) 


-VCft-oO 

(4d) 
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A first step might be the substitution of equation (4e) into equation 
(4b) to give: 

-»[V (ft Fx«+"^0-i) .5, 

Equation (4a) may now be used to eliminate ^ in the above expression 
resulting in: 

DC^-«0“ m'V'Cd- A) ■ •i-Oe’Vwv^Ci- {) (.. 

The above expression may be altered by introducing the following non- 
dimensional force coefficients 

Co- '>(Vv‘s)’‘ 171 

The coefficient may be expanded in a Taylor series as 

C*- Ci. . V q (1^) <•> 

Equation (6) may now be written in terms of and as 

It is now possible to simplify equation (S) somewhat by the following 
redefinitions: 
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(10) 



Equation (10) now allows equation (9) to be rewritten and then rearranged 


as 


Equation (4c), the moment equation, may be written as: 


i 


f‘^%- (frsixHx. 


( 12 ) 


where M has been replaced by C (psiV /2). Again replacing C by a 


Taylor ^ series in a, a 
gives for equation (12) 


q, and 4 using the starred quantities 


0'Jk3‘cA-<^‘c*- ,13, 


Equations (11) and (13) are now the basic equations of planar motion. 
The final goal remains to eliminate one of the variables between theoe 
simultaneous equations. For the present purposes the variaMe q will 
be eliminated and a single differential equation of motion in a will 
be written. As might be expected, this single equation is quite com- 
plicated. Before presenting thlc dynamic equation in a tractable form, 
an outline of the procedure will be given. A fairly straightforward 
approach is to eliminate 4 between equations (11) and (13). The 
resulting equation containing q, a, and a is then differentiated to 
give an expression in and H. Returning to equations (11) and 

(13), eliminating this time q between them now provides a second 
expression q, a, a, and a. Elimination of ^ between these equations 
gives the single dynamic equation in a, a, and a. In carrying out 
the above manipulation it is necessary to perform the different'’ ation 
of (J/V) . This operation may be written as. 




(14) 


The single equation in a that will represent dynamic planar motion 
will be written as: 


where 


H. 


K,’ 








(15) 


(16a) 


(16b) 




(16c) 
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Before attempting to simplify equations (16) it is necessary to reduce 
equation (15) to an equation with constant coefficients as the presence 
of (V/l) introduces time into the coefficients. This may be accomplished 
by writing 

if ■ jiffi ■ ^)^it ■ (i) 

where differentiation has been changed from time, t, to a non-dimen- 
siona'' arc length, (X/i ) . In a similar fashion K may be written in 
terms of (X/l) as 


. (sy- 3KV 


(17b) 


By replacing time derivatives by arc-length derivatives, equation (15) 
now becomes a second order constant co'^fficient equation: 


o< 4- M.« - 1^, 4- R. (if 


(18) 


Admittedly, equations (16) are quite complex; for certain applications 
such as aeroballistics, great simplifications may be made. 

However, before considering this aspect of the problem, the conditions 
for stcd^ility of motion will be examined. 

STABILITY CONSIDERATIONS 


Equation (18) may be rewritten as 


where 


V 2 Xoc' + 

X= 

OlJ = - M, 




(19) 

(20a) 

(20b) 

(20c) 

(20d) 


The term, X, is the deunplng factor of the airship, the term, oiu, is the 
undamped natural frequency. Only for small values of A does the 
body oscillate at this frequency; in the presence of a significant 
amount of damping the planar oscillatory frequency, is less than 
the undamped frequency, The damped planar frequency, u)^, can be 

expressed In terms of A and as 

m - ^Vo/ri) <2i) 

A 

The term «« in equation (20c) is the trim angle of attack due to aero- 
dyneunlc asymmetries, while af in equation (20d) is the trim angle of 
attack due to gravitational path curvature. With regard to the term 
a^, it might be of interest to note that if the airship is neutrally 
buoyant, l.e. s * 1, then from equation (10) § * o and hence from 


// • ■!: 




{■ 


182 



fl . .1 


equation (20d) , ao^ must be zero. 

There are two conditions for oscillatory motion. These are: 

A<.o/w 


(22a) 


> O (22b) 

The former condition allows us to write for the equality A = X = a , 

c n 

(23a) 

where (X/X ) is often called the damping ratio. 

The condition for oscillatory motion given in equation (22b) is that 
0 )^ is real which/ in turn/ from equation (20b) requires that 

< O (23b) 

Under the condition where equations (22b) and (23b) are satisfied/ 
stability (subsident motion) requires that 

;tAs. 

(24a) 


-M. >0 


(24a) 

(24b) 


Thus to assess dynamic and static stability (equations (24a /b) respec- 
tively) it is necessary to assign numerical values to the derivatives 
contained in equations (16a) and (16b) . Numerical values for these 
terms are contained in Table I. While these values may vary with 
airship dimensions/ they have been computed for the airship shown in 
Figure (2) below. An outline of the computational technique is given 
subsequently. 

' Ch.-c^T 

-C* 


Inserting values from Table I gives/ 
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C-2.Z\)C-U^)- C\ + 2 2 07)0 vl <ia) 
-.0437 

-.2 3-1 - 0 437 - -278 


Quite obviously the inequality of equation (24a) is not met so the air- 
ship is not di*namically stable. 

In considering the oscillatory frequency relationship for static 
stability (equation (24b)) we may write: 

It can readily be shown that since the term in the braces is multiplied 
by C* it is rather small. This allows the above expression to be numer- 
ically evaluated as 


M,= l-TS 

< 

I 

Obviously the second condition of equation (24b) is not met. 

It might be expected that nvunerical values of the stjibility derivatives 
would vary from airship to airship. However, it would appear that no 
general simplifications may be made in the and coefficients 
except to omit >erms multiplied by C^. The results seem to indicate 
that for a satisfactory description Jf planar dynamics it is necessary 
to calculate the eight stability derivatives of Table I. Drag, as we 
have noted, is relatively unimportant for estimating the planar dynamics. 


In passing it might be of interest to examine the equivalent expres- 
sions for Hj^, Mx, Ax, and Gx» v’hich are satisfactory for an HTA vehicle. 
If quantities such as Cjg, C^„, and are ignored along with 

the product of starred quantities one, has 
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- vcc 


m. 


(25a) 

(25b) 

(25c) 

(25d) 


Quite clearly the criteria of equations (24) are met when 
and are negative for the HTA vehicle. An examination of ^ 

(16a) and (16b) quickly show that dynamic stedaility cannot depend upon 
such simple criteria in the case of an LTA vehicle: stability conciJ- 
erations are far more complex for the LTA vehicle. 


+ c* • 

equa?? ons 
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For a typical airship we have seen that the motion consists of one 
exponentially undamped mode and one exponentially damped mode since 
from equation (19) 


-X 


and using X = -.139 and Mi = 1.73 we obtain 

X^ « .479 ^2 * ' 

As is well known / the fixed-wing HTA vehicle usually evidences two 
dcunped oscillatory modes. 

CALCULATION OF AERODYNAMIC LOADS 

A fairly straightforward method of calculating static and dynamic loads 
on an airship is the method of virtual mass. While this technique has 
its origin in the work of nineteenth-century hydrodynamicists , it has 
been applied with some success by Bryson to HTA vehicles. Since 
space limitations do not permit even an outline of the derivation, 
reference should be made to either Bryson's work '2) ©r the more readable 
treatment of Nielson '3'. 


Through the use of this virtual mass technique it may be shown that the 
derivatives used in the previous expressions for Hi and Mi are given as, 


where 


C -.ZBCf-ZA, 

C,J-. “1C. 

^ 4 “ _ 

fi, . 


for body-alone and that 




wJ 


(27a) 

(27b) 

(27c) 

(27d) 

(27e) 

(27f) 

(27g) 

(27h) 


(28a) 


(28b) 


for the body in the presence of fins. a(x) is the body radius as a 
function of body station and s (x) is fin span (center-line to tip) as 
a function of body station. In addition Bn, Cn and Dn are defined 


C^im) 


(29a) 
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(29b) 


when "n" and "b" refer to nose and base respectively. 

The above integrals have been evaluated numerically for the airship 
shown in Figure (2) from tabular values of a(x) and s(x). The calcula- 
tions of equations (27) were carried out to give the results shown in 
Table I. 



LENGTH (FT) 

DIAMETER (FT) 

VOLUME (FT3) 

M0M.0FMERTIA(SLU6FT2) 

517 

120 

4.4X10«( 

1.42X10^3 


PIG. 2 REPRESENTATIVE AIRSHIP 


TABLE I 












\d^ 


•730S 

-.o7S6 

+•• 075*6 



Also calculated in the program is the airship volume V = 4.4158 x 10^^. 
Assuming neutral buoyancy it can be shown that (pSL/2m) , equals 1.57, 
which together with values in Table I allows the starred derivatives 
(equation (10)) to be calculated. 

CONCLUSION 

This paper has taken a brief look at the hydrodynamic complexities of 
planar aynamics of airships. It has been shown that the equations of 
motion for a LTA vehicle are far more complex than the corresponding 
, equations of a HTA vehicle. A method has been presented for calcula- 
ting all loads (except drag) acting on a moving airship. 
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